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1. Introduction 

The computation of reliable phenomenological numbers on the basis of 
non-perturbative QCD is still an unsolved problem. Numerical methods 
seem to be required, as the simulation of the theory on a space-time lattice. 
A semi-phenomenological tool is provided by the Schwinger-Dyson equations 
(SDEs), which are reviewed in this respect in [1]. In the present paper we 
discuss an approach to non-perturbative phenomena in non-abelian gauge 
theories, which is based on the integration of exact renormalization group 
equations (ERGEs) [2] in continuum quantum field theory [3]. 

ERGEs describe the continuous evolution of effective Lagrangians, or ef- 
fective actions, with a scale (or infrared cutoff) k, and allow to obtain the 
full quantum effective action for k — > from a bare (microscopic) effective 
action at a cutoff scale k = A. In contrast to standard renormalization group 
equations they describe this evolution including all irrelevant couplings, or 
higher dimensional operators, and are exact in spite of the appearance of 
only one-loop diagrams. The price to pay is the fact that one has to deal 
with an infinite system of coupled differential equations which for practical 
purposes requires some truncation. The types of approximations or system- 
atic expansions, which can be employed, depend on the phenomena under 
consideration: 

One can expand the effective Lagrangian (or Hamiltonian or action) in 
powers of momenta or derivatives, but keep all powers of the involved field. 
This kind of expansion within the ERGEs is appropriate for the calculation 
of effective potentials, which are required to study the existence and nature 
of phase transitions (see [4] for some early literature). 

Alternatively, one can expand the effective action in powers of fields, 
keeping all powers of the momenta [5]. This allows, e.g., to study the forma- 
tion of bound states by using ERGEs to look for the appearance of poles in 
four-point functions [6]. Finally, in order to describe dynamical symmetry 
breaking, composite fields can be introduced [7, 8] and different expansions 
can be applied to the parts of the effective action involving fundamental and 
composite fields, respectively. 

The ansatz of the present paper is based on the expansion of the effec- 
tive action in powers of fields, which allows to study the full momentum 
dependence of the (gluon and ghost) propagators. Within this expansion the 
knowledge of the running 3- and 4-point functions is required in order to 
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integrate the ERGE for the running of the 2-point function, the knowledge 
of the 4- and 5-point functions is required for the ERGE of the 3-point 
function etc. [5, 6]. 

In the case of gauge theories this expansion gets to some extent reorga- 
nized: due to gauge invariance (or, more precisely, BRST invariance) parts of 
higher N-point functions (with certain Lorentz and gauge index structures) 
are completely determined by lower N-point functions. The usual procedure 
consists in imposing Slavnov-Taylor identities (STIs) on the effective action, 
which can be solved for "dependent" parts of N-point functions in terms of 
"independent" ones. For instance, in the present case of a pure Yang-Mills 
theory a gauge invariant term in the effective action of the form 

f;m{d 2 )f, u t (i.i) 

(where D denotes the covariant derivative and / an arbitrary function) de- 
scribes simultaneously the "independent" effective gluon propagator or 2- 
point function, and certain "dependent" Lorentz and gauge index structures 
of the 3-, 4- and higher gluonic N-point functions. Additional "independent" 
structures of the gluonic 3- and higher N-point functions are described in 
terms of gauge invariant expressions involving at least three powers of the 
field strength F M „. 

The application of ERGEs to gauge theories faces some technical prob- 
lems, since ERGEs are based on an intermediate infrared momentum cutoff 
k, and such a cutoff generally breaks gauge or BRST invariance [9-14]. One 
way to deal with these problems is the modification of the STIs [11, 13, 14]: 
the modified STIs impose "fine tuning conditions" on those couplings in the 
effective action at k ^ 0, which break gauge or BRST invariance. These "fine 
tuning conditions" are such that they guarantee the BRST invariance of the 
effective action for k — > 0. 

As a consequence of the (modified) STIs, at any value of k, the "de- 
pendent" parts of the effective action can thus be obtained either via the 
modified STIs in terms of the "independent" parts at the scale k, or directly 
from the integration of the ERGEs. Below we will employ approximations 
which consist in neglecting contributions to the r.h. sides of the ERGEs. In 
the absence of approximations, both methods to obtain the dependent parts 
lead to identical results. In the presence of approximations one should de- 
termine the dependent parts of the effective action from the modified STIs. 
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Approximations on the r.h. sides of ERGEs for independent parts of the ef- 
fective action do thus not imply a violation of the modified STIs. 

In the presence of approximations one can nevertheless study the ERGEs 
for the dependent parts of the effective action and see, whether these ERGEs 
deviate strongly from the corresponding equations obtained from the k- 
derivatives of the modified STIs. As emphasized in [14], this method allows 
for a numerical self-consistency check of the approximation and will also be 
employed in the present paper. 

As described in section 3, our approximation consists in the neglect of the 
contribution of certain parts of the 3- and 4-point functions to the ERGEs 
of the 2-point functions. We take those parts of the 3- and 4-point functions 
into account, however, which are determined in terms of the 2-point functions 
through the standard part of the STIs. In the purely gluonic sector these are 
those obtained by an expansion of eq. (1.1) in powers of the gauge field A a . 

No approximation is imposed on the momentum dependence of the 2- 
point functions (the gluon and ghost propagators). These propagators by 
themselves are gauge dependent and not directly observable quantities. A 
certain combination thereof, however, has a physical meaning: The potential 
between infinitely heavy quarks can be written as a product of the dressed 
quark-gluon vertices and the full gluon propagator, and the quark-gluon ver- 
tex function is given by the ghost propagator function. As a result the heavy 
quark potential depends on both propagator functions, which we calculate 
here in the quenched approximation. 

The aim of this paper is to see, how much information on the heavy 
quark potential in quenched QCD can be obtained within the approximation 
we employ, which is seen to be self-consistent for a range of momenta p 2 
down to a certain value p 2 ~ /c 2 (see section 5). We find strong evidence for 
a l/p 4 -behaviour within the trustworthy range of p 2 . 

Both the method and the results can be compared with the Schwinger- 
Dyson formalism. To this end we derive in section 2 a new and general 
relation between ERGEs and SDEs: we show that effective actions, which 
satisfy SDEs in the presence of an infrared cutoff k, are (quasi-) fixed points 
of the ERGEs. Thus a search for fixed points of the ERGEs for k — > is 
equivalent to the search for solutions of SDEs. 

Investigations of the gluon propagator have been performed in the con- 
text of SDEs both in the axial [15, 16] and Landau [17] gauge. In most 
cases solutions with a singular behaviour like l/p A have been found to be 
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consistent. (A behaviour less singular than 1/p 2 has been put forward on 
phenomenological grounds in [18] and claimed to exist in the axial gauge in 
[16], but recently such a form of the gluon propagator in the infrared has 
been argued to be inconsistent [19].) In the Landau gauge the ghost sec- 
tor has been completely neglected; thus it is much less straightforward to 
relate these results to physical observables. In section 6 we will discuss the 
comparison between the ERGE and SDE approach in more detail. 

The paper is organized as follows: in section 2 we discuss the derivation of 
the ERGEs, and their relation to SDEs. In section 3 we present our truncated 
ansatz for the pure Yang-Mills action. In section 4 we discuss the technical 
procedure, and in section 5 our numerical results. Section 6 is devoted to 
conclusions and an outlook. 

2. ERGEs and SDEs 

The form of the ERGEs can most conveniently be derived assuming a 
path integral representation for the generating functional of the theory. To 
this end we restrict ourselves, for the time being, to the case of a single scalar 
field. It is straightforward to supplement the following equations with the 
appropriate minus signs, Lorentz and internal indices in the cases of fermions, 
vector fields and internal symmetries. To start with, the path integral for- 
mulation of the theory requires for its proper definition a regularization of 
the ultraviolet divergences. This will be implemented through a modifica- 
tion of the propagator, i.e. the term quadratic in the fields in the bare action. 
In addition we introduce an infrared cutoff k; the ERGEs will describe the 
variation of the generating functional with k, where the UV cutoff A will be 
held fixed. (In the framework of the ERGEs, in particular in the application 
to gauge theories, one can consistently set A = oo, see below). The UV and 
IR regularized generating functional of connected Green functions Q^{J) can 
then be represented as 



where Si nt (<p) is independent of k. Here (J, ip) etc. is a short-hand notation 




{ V ,P^ V )-AS^-S mt { V )+{J,v) 



(2.1) 



for 




(2.2) 
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P 1 denotes the bare inverse propagator (e.g. P 1 {p 2 ) = p 2 in the case of a 
massless scalar field), S int the bare interaction, and 

U3t = \{tp,Bbtp) (2.3) 

implements the UV and IR cutoffs. In the case of a massless field a convenient 
choice for R% is 

2 2 



e a- 



such that the full propagator 



J? = (P.- 1 + *£)" = ' - / ' (2.5) 

has the desired properties: exponential decay in the UV, and finiteness in 
the IR. 

The ERGE for Qfc describes its variation with k, at fixed A. It can be 
obtained by differentiating eq. (2.1) with respect to k, and replacing ip under 
the path integral by variations with respect to the sources. One finds 



\ s 2 g£(j) sgt(J) sgt(J) 

SJ(q) SJ(-q) 5J(q) SJ(-q) 



(2.6) 

The effective action r£(ip) is defined by the Legendre transform of Q^(J), 

T£(V)=G£(J) + UV) , (2.7) 

and it is convenient to define an effective action with the IR cutoff terms 
subtracted: 

f £ = It - AS fc A (2.8) 

Inserting the Legendre transformation into eq. (2.6) one obtains the ERGE 
for f£: 
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The boundary condition of 1^(93) for k — > A can be obtained from a careful 
consideration of the limit — » for /c — > A and reads 

Tl( V )= l -( ! p,{P£r 1 V )+SUv) • (2.10) 

On the other hand, for k — > r£(</?) becomes the full quantum effective 
action of a theory, which is defined by a bare interaction S int (<p) and a fixed 
UV cutoff A. Since the integration of the ERGEs allows us to relate effective 
actions with different k, it can be used to compute Tq in terms of 
related to Sm by (2.10). 

For most applications in particle physics this is actually all we want: 
there always exists a scale A (at least the Planck scale) beyond which we 
expect unknown new interactions or particles with masses larger than A. 
Physics beyond A, including quantum fluctuations involving momenta p 2 
with p 2 > A 2 , will certainly affect the form of Si nt or the effective action T\. 
Given an ansatz for 1^, however, we are generally interested in the effect of 
quantum fluctuations involving momenta with < p 2 < A 2 . If we integrate 
the ERGE for T% from k = A down to k — 0, these quantum fluctuations 
have been fully taken into account. 

Next we turn to the relation between ERGEs and SDEs. To this end it 
is first convenient to write the ERGEs in terms of the generating functional 
(or partition function) Zfc(J), where Z£(J) is given by 

Z fc A (J) = e- g ^ (2.11) 

and the r.h.side of eq. (2.11) by the path integral (2.1). The ERGE for Z£(J) 
reads 

^"l/^^jTO ' (2 - 12) 

The full set of SDEs can also most conveniently be expressed in terms of 
Z(J) [20]. They follow from the vanishing of the path integral over a total 
derivative with respect to tp. Subsequently we will consider the SDEs in the 
presence of an UV cutoff A and an IR cutoff k as in (2.1). In terms of Z£(J) 
they read 



Ap)-(Po\p 2 ) + RUp 2 )) 



SJ(-p) 5tp(-p) 



V = Ti, 



z£(j) = o . 

(2.13) 
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After replacing Z£ by Qfc according to eq. (2.11), and Q£ by according 
to eq. (2.7), one obtains the SDEs for the one-particle irreducible vertices, 
which are generally considered in the literature for k = and A = oo. 

Let us now denote the entire l.h.side of eq. (2.13) by E£(J). Next we 
evaluate the fc-derivative of E£. The /c-derivative hits R% and Z£, and we 
have to use the ERGE (2.12). The result can be written in the form 

Thus we can make the following important statement: suppose that we man- 
aged to find a functional Z£(J), which satisfies the SDEs (2.13) for some 
value of k, identically in the sources J. Thus the corresponding quantity 
Efc(J) vanishes. Now let us start to vary k. Since the r.h.side of eq. (2.14) 
vanishes, Ejjh will still vanish for k' ^ k. (Here, of course, it has been used 
that Z£ varies with k according to the ERGE (2.12).) Thus we find that 
Zu still satisfies the SDE (2.13), with k replaced by k' . In other words, if 
we found a functional Z k which satisfies the SDE (2.13) for some value of 
k, and compute Zfe by integrating the ERGE (2.12), Z§ is guaranteed to 
satisfy the SDEs for all values for k'. The same statement will hold for the 
other generating functionals Q£ and r£. Thus generating functional which 
satisfy /c-dependent SDEs can be named quasi-fixed points of the ERGEs 
("quasi" because they still depend on k). In particular, for k — > 0, Qq and 
Tq will satisfy the standard SDEs (without an infrared cutoff, but still an 
ultraviolet cutoff A) if such a quasi-fixed point is approached in the infrared. 

Note that this property is independent of the form of S int in eq. (2.13): 
Any SDE with any form of Si n t in eq. (2.13) plays the role of a quasi-fixed 
point. (In particular, since we work in the presence of an UV cutoff A, Si n t 
can contain arbitrary higher dimensional interactions suppressed by powers 
of A.) 

These formal arguments let us suspect, that results obtained by the in- 
tegration of the ERGEs can be quite similar to results obtained by solving 
SDEs, although the formal arguments no longer hold in a strict sense as soon 
as the effective actions are approximated or truncated in some way. 

Let us note, for later use, that on the r.h. sides of the ERGEs (2.6) and 
(2.9) we can replace R k (p 2 ) by R k (p 2 ), where R k (p 2 ) is defined by 

R k (p 2 ) = Rt°°(p 2 ) • (2-15) 
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In the limit A — > oo the momentum integrations on the r.h. sides of (2.6) and 
(2.9) are still exponentially damped with the choice (2.4) for R k . Accord- 
ingly we denote by Q k and Y k the generating functional, which satisfy the 
ERGEs with replaced by R k . Of course we are still able to compute T k=0 
in terms of an ansatz for T k= \. In physical terms UV divergences do not 
appear, since I\ is already supposed to contain all quantum effects involving 
momenta with p 2 > A 2 . We just loose the possibility to formally equate I\ 
with Si n t appearing in a path integral of the form of eq. (2.1), although a 
slightly different, more involved path integral representation for T k still ex- 
ists [21]. From universality we expect in any case, that Tq is independent 
of minor modifications of F\, so one choice is a priori as good as the other. 
Thus, subsequently, we define a theory by the ERGEs (2.6) or (2.9), with R k 
replaced by R k , and a finite ansatz for I\ (or Q\) at a fixed scale A. 

3. Yang-Mills Theories 

The application of ERGEs to gauge theories requires some treatment of 
the breaking of gauge invariance, which is introduced through the interme- 
diate IR cutoff k. Here we use the method of imposing modified STIs [11, 
13, 14] in the formulation presented in [13, 14]. In ref. [14], in particular, 
ERGEs for SU(3) Yang-Mills theory have already been integrated numeri- 
cally within an approximation, where only marginal and relevant couplings 
in the effective action have been taken into account. 

In the present paper we extend the form of the effective action consider- 
ably, namely we allow for an arbitrary momentum dependence of the gluon 
and ghost propagators. Again we consider the SU(3) Yang-Mills theory in 
four-dimensional Euclidean space-time, where the classical action reads 

S = J & A x [\F; V F; U + ^p.A^Al + d,c a (d,c a + gf\ c Ay) 

-K; (d,c a + gf\ c Ay) - L a \gr bc c h c c + L a U,A a ^ (3.1) 

with 

F; v = d,A a v - d v Al + gj\cA\Al . (3.2) 

We have included the usual gauge fixing and ghost parts and coupled external 
sources to the BRST variations 
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5AI = (d,c a + gf\cA\c c ) ( 
5c a = \gf\ c c b cX 



Sc a 



a 



(3.3) 



where ( is a Grassmann parameter. The invariance of 5 at L = under 
BRST transformations can then be expressed as 



y ( M 5A« 5c a 5c a J 

= c /dUi^ 



55 55 55 55 



5L a 5c a 5L a 5c a 



L=0 



L=0 



(3.4) 



IR cutoffs in the gluon and ghost propagators will be introduced by adding 
a term of the form 



1 



AS k = -(A«,R k ^ A a v ) + (c\R k c a ) 



(3.5) 



to the action. Explicit expressions for R k:flI/ and R k will be given below. The 
effective action 

r k (A,c,c,K,L,L) =f k + AS k (3.6) 
satisfies an ERGE of the form [13, 14] 



^=/(S{^^ 2 )(^):: ( . pMS 

- d k R k ( P 2 ) (r 



(p) 



,(2) 
k 



-1 

-c a (-p),c a (p) 



(3.7) 



Here (r^)^ ^ Ab( . ^ is the p), A* (g))-component of the inverse of the 



matrix rjj. 2 ^ = 5 2 T k /5(p 5(p of second derivatives of T k with respect to the fields 
(p B = {A®(—p),—c a (—p),c a (—p)} and tp B = {A a {p), c a (p), c a (p)}, where the 
index of <p and tp runs over the different fields, momenta, and Lorentz and 
gauge indices. For k — > 0, T has to satisfy the STI (3.4). To this end T k has 
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to satisfy, in particular at the starting point k = A, a modified STI of the 
form [13, 14] 



d 4 p \ ST k 6f k 6f k 5f k Sf k Sf k 

(2tt) 4 \6K£(-p) 5A a ^(p) 5L a (-p) 5c a (p) 5L a (-p) 5c a (p) 



L=0 



~ J (2tt)4 ^ [ Kk ^ [P } 8K-{-p) 5 V B ^ k ) 

~ / 2\ 5 2 Y k / (2)\ — i 

~ Rk{jP >5L a {-p) 5<p B v Ffc h B ,c«ip) 

The crucial point is that once eq. (3.8) is satisfied for k — A, it will be 
satisfied by T k for any k < A provided Y k is obtained from T\ by integrating 
the ERGE. In particular f fc=0 will satisfy the standard STI (3.4), if the IR 
cutoff functions R k ^ u and R k vanish identically for k — > 0. 

Most importantly, T k at the starting point k = A can thus not be identi- 
fied with the classical action (3.1), but it has to contain symmetry breaking 
terms specified by the need to satisfy eq. (3.8). The origin of these sym- 
metry breaking terms can be understood from the interpretation of I\ as a 
quantum effective action, where momenta p 2 with p 2 > A 2 have already been 
integrated out; the symmetry breaking induced by this "IR" cutoff A then 
generates the symmetry breaking terms in T\. 

As stated above, our approach towards the integration of the ERGE (3.7) 
is based on an expansion in powers of the gluon and ghost fields. We are 
interested in the form of the gluon and ghost propagators, since a certain 
combination thereof (see below) determines the heavy quark potential in 
the quenched approximation. It is thus essential to allow for an arbitrary 
momentum dependence of the 2-point functions in our ansatz for the effective 
action. 

Expanding the ERGE (3.7) in powers of fields one easily finds that the full 
3-gluon, 4-gluon, ghost-gluon and 4-ghost vertices appear on the r.h. sides 
of the ERGEs for the gluon and ghost propagators. In the sequel we will 
refer to parts of the vertices with certain Lorentz and gauge index structures 
as "operators" . One can distinguish two different kinds of such operators: 
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a) "Independent" operators, which are not fixed in terms of the prop- 
agators or other operators by the STL An example is the contribu- 
tion to the 4-gluon vertex obtained from a term involving four powers 
of the field strength F£ v in the effective action. (The r.h.side of the 
corresponding ERGE depends on the full 2-, 3-, 4-, 5- and 6-point 
functions.) 

b) "Dependent" operators, which are fixed by the modified STI in terms 
of the propagators and operators of type a). In general eq. (3.8) yields 
a system of coupled non-linear integral equations for the dependent 
operators. Examples are the operators corresponding to those parts of 
the 3- and 4-gluon vertices, which are obtained by an expansion of the 
expression F^ u {f{D 2 )F^ u ) a in powers of the gauge field (see eq. (1.1)), 
and a gluonic mass term (which we do not consider as part of the gluon 
propagator in this context). 

Since our goal is the determination of the propagators, the only impor- 
tant operators for our purpose are those, that contribute significantly to the 
ERGEs for the 2-point functions. In a first approximation, we will neglect 
the contributions of the operators of type a), which are the ones not de- 
termined by the STI. Such an approximation is only expected to hold for 
a limited range of scales k, and can to some extent be checked a posteriori 
(see below). In any case, it turns the ERGE (3.7) into a closed system of 
integro-differential equations for the propagators, which can be integrated 
numerically. 

Furthermore, with the simpler ansatz used in ref. [14] it has proved to 
be permissible in a quantitative sense to determine the operators of type b) 
through the standard form of the STI (eq. (3.8) with a vanishing r.h.side) 
instead of using the modified STI, with the important exception of the gluon 
mass term. Again, this may only be true within a certain range of k. We 
will therefore proceed in the same way here, i.e. we assume that we can 
approximate the operators of type b) by the corresponding ones obtained 
from the standard STI, as far as their quantitative contribution to the flow 
of the 2-point functions is concerned. 

In practical terms, this assumption implies a tremendous technical sim- 
plification in the determination of these operators. For example, for the 3- 
and 4-gluon vertices we can use the expressions resulting from an expansion 
of F® u (f(D 2 )F[ lu ) a in powers of A® (apart from a momentum-dependent field 
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renormalization, cf. the explicit expressions in the appendix). We consider it 
to be essential, however, to use the modified form of the STI to fix the gluon 
mass term, which is a truly relevant parameter of the effective action. 

Let us stress that the approximation we employ for the operators of type 
b) does not imply that we are working with an effective action, which violates 
the modified STI. In principle we could, for every value of k ^ 0, determine 
the operators of type b) by solving the system of coupled non-linear integral 
equations given by eq. (3.8), so that the resulting effective action would be 
in accord with the modified STI. What we actually do is to approximate the 
contributions of the operators of type b) to the flow of the 2-point functions 
by the contributions of the corresponding operators obtained from the stan- 
dard STI, in the same way as we neglect the contributions of the operators 
of type a). 

It is important to note, that there exists a non-trivial quantitative check 
for the self-consistency of the approximations employed. It is based on the 
consistency of the ERGE with the /c-derivative of the modified STI in the 
absence of approximations, and will be discussed in more detail in the next 
section. In particular, this self-consisteny condition will be used to determine 
the range of scales k, where we consider the approximations to be permissible. 

The approximations employed here can be improved systematically by 
taking more and more operators of type a) into account, as well as by using 
the operators of type b) corresponding to a (partial or complete) solution 
of eq. (3.8). The purpose of the present paper is to see, which results for 
the heavy quark potential will be obtained within the "minimal" approach 
described above. 

Let us now write down those terms in the /c-dependent effective action, 
which we take into account on the r.h. sides of the ERGEs for the 2-point 
functions. In a slight abuse of notation, we will denote this part of the effec- 
tive action by T k . Firstly, we include the gluon and ghost 2-point functions 
themselves, which are given by the momentum-dependent functions and 
fi,ki respectively, and where we take a mass term for the gluons into account. 
Since we will consider the Landau gauge later on, it is not necessary to in- 
troduce a further momentum-dependent function in the longitudinal part of 
the gluon 2-point function here. Secondly, we describe the 3-gluon, 4-gluon 
and ghost-gluon vertices of the type b) above (determined from the stan- 
dard STI), which we take into account on the r.h. sides of the ERGEs for the 
2-point functions. Finally we include the external source terms containing 
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K®, L a and L a , which are required to formulate the STI. These latter terms 
are constrained by two further identities satisfied by T k , which can be shown 
to be invariant under the RG flow [11, 13, 14]: 

§ = ^Al . (3.9b, 
Explicitly, then, our ansatz for T k reads 

r k (A,c,c,K,L,L) = 

\ ! a«(pO { ( P j v - Pl » Plv ) f\M) + + ™* <wj ^(ft) 

+ i £fc Ac / /•WPi,P2,P 3 ) ^(Pi) ^(pa) ^(p 3 ) 

Jpi,P2,P3 

+ £fc Ac Ae / U„(Pl,P2,P3,P4) A*(p 2 ) ^(p 3 ) ^(p 4 ) 

Jpi,...,p 4 

+ / c a ( Pl ) p? / 2>fc (p?) c a (p 2 ) 

Jpi,P2 

+ < <?* Ac / C a ( P l) Pi,, / 2> *(P?) ^(rf) 4(P2) C C (p 3 ) 

+ * / p lM / 2>fc (p?) c a (p 2 ) 

Jpl,p2 

- g k f\ c I k;( Pi ) f 2 M) KM) Ab M) cC fe) 

Jpi,P2,P3 

- 9 i Ac I L\ Pl ) c\p 2 ) c c (p 3 ) 

^ J Pl,P2,P3 

- % -\ L a ( Pl ) Pl , A;(p 2 ) , (3.10) 
a Jpi ,p 2 



>Pl,P2 

where 
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The /c-dependent 3-gluon and 4-gluon vertex functions fz^ vp and f±^ vpiJ are 
given in terms of the gluon propagator function /i^ and the ghost propagator 
function / 2 fe in the appendix. 

In order to establish the connection with the heavy quark potential, we in- 
troduce quark fields in the effective action. In the limit where the quark mass 
tends to infinity, all quantum corrections corresponding to diagrams which 
contain inner quark lines, are naively suppressed by powers of the quark mass. 
On the other hand, the STI requires a momentum dependent quark-gluon 
vertex function, where the momentum dependence involves the function ji,k 
above. (Within the standard approach, the origin of this non-trivial contri- 
bution is due to UV divergent diagrams, which are not suppressed by the 
heavy quark mass.) After an eventual redefinition of the quark fields, such 
that their kinetic term is properly normalized, the standard STI constrains 
the quark-gluon coupling in the effective action to be of the form 

- 9k \ a AB I ¥(Pi) 1, UM) A«(p2) ¥(Ps) , (3.12) 

J Pl,P2,P3 

where A denotes the SU(3) generators in the fundamental representation. 
Here we have implicitely used the i^c-coupling from eq. (3.10) and the fact 
that the quantum corrections to the vertices, which couple external sources 
to the BRST- variations of the quark fields, are suppressed by powers of the 
quark mass in the Landau gauge. 

Let us now consider the 2-point function of the quark currents J°, 
given by 

J a , = *ab ¥ 7, r • (3.13) 

In the heavy quark limit, the only contribution to comes from the dressed 
one-gluon exchange diagram, so from eq. (3.12) and the form of the gluon 
propagator we obtain the current 2-point function 

ytf) ^ww ' (3 - 14) 

We have omitted the gluon mass term m| in our definition of Vt(p 2 ). In 
principle we could have replaced the factor p 2 /i,fc(p 2 ) in the denominator of 
^4(p 2 )> which is due to the gluon propagator, by the factor involving to 2 , 
and the infrared cutoff term. However, once the integration of the ERGE 
is pursued down to k — 0, the modified STI (3.8) turns into the standard 
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STI and enforces m| =0 = 0, and the infrared cutoff term vanishes identically. 
Using eq. (3.14) to define Vk(p 2 ), these effects have been anticipated, and 
V k (p 2 ) corresponds more closely to the physical function V (p 2 ) already at 
finite k. 

In the heavy quark limit V k (p 2 ) can be identified, up to trivial factors, 
with the Fourier transform of the potential between quarks in quenched QCD. 
Although we will compute the gluon and ghost propagator functions f i;k 
and /2,/fc individually by integrating the ERGEs, it is only the combination 
appearing in eq. (3.14), which has a physical meaning. 

A priori it may seem that the functions fi yk and f2, k can be changed at 
will by redefining the gluon, ghost and anti-ghost fields in the form of multi- 
plication with momentum dependent functions. We have checked explicitely, 
however, that the STI restricts the relations between the vertices and the 
two-point functions such that finally all possible field redefinitions cancel in 
the expression for V k (p 2 ). 

As noted above and in refs. [11, 13, 14], ml in is fixed in terms of the 
other parameters in Tk by the non-vanishing r.h.side of eq. (3.8), where we 
use our ansatz (3.10) for T k . After a consideration of terms ~ A ■ c of eq. 
(3.8), one arrives at a non-linear relation of the form 

m 2 k = j^- 2 k 2 ST(k 2 ,m 2 k ,f 1>k ,f 2 , k ) , (3.15) 

where we have anticipated the Landau gauge and omitted the gauge param- 
eter a on the r.h.side. The expression ST(/c 2 ,...) involves momentum inte- 
grations over a number of one-loop diagrams, in which the gluon and ghost 
propagators and hence the functions fi yk and appear, and will be given 
in the appendix. Generally eq. (3.15) can be solved for mf only numerically, 
but in any case m| is not an independent parameter of Tfc. 

To determine the flow of the coupling parameter g k , we will use the ERGE 
for the KcA- vertex with vanishing momenta. Then it turns out, that in the 
Landau gauge g k actually does not run at all with k; it remains a constant, 
which will be specified once and for all at the starting point k = A. (It is 
straightforward, on the other hand, to define a "physical" running coupling 
in terms of the dressed one-gluon exchange diagram or the function V k of 
eq. (3.14), which runs with k due to the running of fi >k and /^.fc- Within 
our approximation this running agrees to one-loop order with the usual per- 
turbative running, as we have checked explicitely.) In the Landau gauge the 
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independent running parameters in Tk are thus only the two functions f\ t k 
and /2,fe. The aim will be the determination of these quantities by integrating 
the ERGEs with suitable boundary conditions. 

Finally we have to introduce and specify the IR cutoff term ASk in eq. 
(3.5). If we demand a reasonable behaviour of the full propagators including 
the IR cutoff terms which appear on the r.h.side of the ERGE (3.7), we are 
lead to more sophisticated choices than an expression of the form of eq. (2.4) 
with A — > oo. (Note that Rk^ v and Rk have to vanish identically for k — > 
in order that the modified STI (3.8) turns into the standard one (3.4). To 
this end A has to be put equal to oo inside Rk^v and Rk, which is perfectly 
consistent as discussed in sect. 2.) First, we do not want to generate any poles 
on the real axis through the introduction of ASk- This enforces a dependence 
of Rk,fiu and Rk on the functions /i^ and / 2 ,fc such that they are proportional 
to fi t k and / 2) fc, respectively. Now, however, it is no longer automatically 
guaranteed, that ASk serves as an IR cutoff if, e.g., fi t k(p 2 ) vanishes for 
p 2 — > (as we might possibly expect). This problem is overcome by writing 
the functions /i^ and $2,k also in the exponents. Taking the presence of the 
gluon mass term m\ into account, our choice of the cutoff functions thus 
reads 



e -(p 2 fiAp 2 )+™ 2 k )/k 2 
RkMP ) = [P hkip )+m k ) - _ e - (p2/l , fc(p2)+mD/fc2 



VvVv 



/p2 \ e -{p 2 /a+m 2 k )/k 2 

+ {a +ml ) l- e-W/^nl)/^ ^ ' ( 3 - 16a ) 

e -P 2 / 2 , fc (P 2 )/fc 2 

Rk(p 2 ) = P 2 f2, k (p 2 ) i _ e - P i hAp i)/k 2 ■ ( 3 - 16b ) 

In order to obtain the full propagators, these terms have to be added to 
the quadratic terms ~ A ■ A or ~ c • c in the action Tk- Let us at this point 
remark a subtlety which has already been noted in ref. [14]: the modified 
STI (3.8) determines, after its expansion to 0(A ■ c), only the longitudinal 
part of the gluon propagator and requires the presence of a term of the form 
m k PtiPv/p 2 - Then it is the condition of locality of the effective action, i.e. 
the need to cancel terms of the form p^Pvjp 2 , which requires the presence of 
the same mass term m| in the transverse part of the gluon propagator. Only 
after this consideration has been completed, we are allowed to approach the 



16 



Landau gauge a — 0. In ref. [14] we had checked that a = is preserved 
by the RG flow (and IR stable), and that the r.h. sides of the ERGEs remain 
finite. Since this choice simplifies the calculations considerably, we will work 
in the Landau gauge throughout the rest of this paper. The full gluon and 
ghost propagators, as derived from T k = T k + AS k , then read 

' V & v 1 (3.17a) 



p 2 fiAp 2 ) + m l \ ^ p 2 

1 _ e -p 2 /2,fe(p 2 )/fc 2 
P 2 f2,k(P 2 ) 



(3.17b) 



and 

respectively. 
4. Procedure 

Let us now turn to the computation of fi >k and f'2,k by integrating the 
ERGEs. The ERGE for fi >k can most easily be obtained by studying the 
terms quadratic in A of the functional ERGE (3.7). Its diagrammatic form 
is shown in fig. 1. Note that the 3- and 4-gluon vertices are more complicated 
than the ones of a classical Yang-Mills action and are given in the appendix. 
The form of the 3-gluon vertex coincides with the one used in the framework 
of SDEs [17] because it corresponds to a special solution of the standard STIs 
[22] (up to a momentum-dependent field redefinition); we have not been able 
to find the corresponding 4-gluon vertex elsewhere in the literature. The 
ERGE for f 2tk is obtained from the terms ~ cc of eq. (3.7), and is also shown 
in diagrammatic form in fig. 1. 

The crossed circles in fig. 1 denote insertions of d k R k ^ u resp. d k R k accord- 
ing to the ERGE (3.7). Due to the complicated expressions (3.16) for these 
cutoff functions, i.e. their dependence on m k , f\ yk and f2,k, the /c-derivatives 
d k acting on R k ^ v and R k have to be decomposed into partial derivatives 
d/dk, df 1>k /dk ■ d/dfi }k , etc. 

The system of ERGEs for f 1>k and f 2 , k (in the Landau gauge), which is 
obtained from the diagrams of fig. 1, is thus of the form 

dkfi,k = (hi(k 2 ,m 2 k , /i,fc,/ 2 ,fc) + d k f hk * h\(k 2 ,m 2 k J 1M f 2 , k ) 

+ d k f 2 , k * hl(k 2 , m 2 k , f hk , / 2)fe ) + d k m\ ■ h\(k 2 , m 2 k , / ljfc , f 2 , k )) , (4.1a) 
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d k f2,k = (h° 2 (k 2 , m 2 k , f 1>k , f 2tk ) + d k f ltk * h\(k 2 \ m\, f 1>k , f 2 , k ) 

+ d k f2,k*hl(k 2 ,mlJ ltk J 2 ,k) + d k ml-hl( y k 2 ,ml,f ltk ,f2,k)) • (4.1b) 

The evaluation of the expressions h\ requires a numerical computation 
of the one-loop integrals appearing in fig. 1. The * denotes corresponding 
convolutions with respect to the loop momentum. It is evident from eqs. 
(4.1), that we still need the knowledge of d k m k in order to determine d k fi^ k 
and d k f 2yk . d k m k can be obtained from the same diagrams in fig. 1, which 
are relevant for d k f\ yk , in the limit of vanishing external momentum. Again 
this gives us an expression of the form 

2 

d k m 2 k = (h° 3 (k 2 , m 2 k , f 1:k , f 2 , k ) + d k f hk * h\(k 2 , m 2 k , f ljk , f 2 , k ) 

+ d k f 2jk *h 2 3 (k 2 ,m 2 k J hk J 2jk ) + d k m 2 k -hl(k 2 ,m 2 k ,f hk ,f 2<k )) . (4.2) 

Alternatively, we can obtain a similar expression by differentiating eq. 
(3.15) with respect to k. We thus have two different equations at our disposal 
to complete the system (4.1). In the absence of approximations, the two 
different completed systems would yield exactly the same results, as follows 
from the compatibility of the full equations (3.7) and (3.8). Since the r.h. sides 
of the eqs. (4.2) and (3.15) involve our ansatz (3.10) for T k rather than the 
full effective action, however, we obtain different results. (Nevertheless, they 
agree on the one-loop level [14].) The use of T k instead of the full effective 
action involves two different kinds of approximations related to the operators 
of type a) and b), respectively (see the previous section). 

The compatibility of the two completed systems of equations considered 
above turns out to be sensitive to both kinds of approximations. As in ref. 
[14], we will use the difference between the corresponding results as a quan- 
titative measure of the inadequacy of the approximations (for the respective 
scale k). Here we concentrate on the effect of the approximations on the 
gluon and ghost propagator functions: We will take a certain discrepancy (in 
fact 10%) in the results for d k f\ tk or d k ] 2)k at some value of the momentum 
as an indication that the approximations have a considerable influence on 
the contributions to the flow of f± tk or f 2>k . 
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We are of course aware of the fact that the described consistency check 
merely yields a necessary condition for the validity of the approximation, 
not a sufficient one. There exists, however, another a-posteriori justification 
for our procedure: The results we obtain will show a remarkable degree of 
universality in the sense of ref. [23] (see the next section) as long as k > k , 
where the scale ko is determined from the above consistency condition. Since 
universality is only expected to hold when gauge symmetry is taken into 
account properly, and no significantly contributing operators are neglected, 
this property provides a strong indication of the consistency of our procedure. 

In practice we will use the form of eq. (4.2) as derived from the ERGE 
(3.7) in order to complete the system of linear equations for dkfi,k and dkf2,k- 
The other system of equations, with eq. (4.2) taken from the derivative of 
eq. (3.15), will be employed for the consistency check at every scale k. 

Given the ERGEs for and /^fc, we still need suitable boundary con- 
ditions in order to proceed with the numerical integration. These boundary 
conditions can be chosen such that T\ resembles as close as possible the 
classical (bare) Yang-Mills action (3.1). One would thus choose, at k — A, 
/i,a(p 2 ) = /2,a(p 2 ) = 1 and not too large in order to start in the pertur- 
bative regime. It is possible, however, to improve the form of the effective 
action at the starting point k = A. By definition, T\ should include all quan- 
tum effects involving internal momenta p 2 > A 2 . On the one hand these are 
certainly small for a small coupling g\, on the other hand it is fairly straight- 
forward to include them to one-loop order. To this end one has to calculate 
the one-loop diagrams, which contribute to the gluon and ghost propaga- 
tor, with IR cutoff functions Rk as in eqs. (3.16) and k = A. (Inside these 
cutoff functions we can set, to one-loop order, fi,k(p 2 ) = f2,k(p 2 ) = 1 and 
m\ = 0). The diagrams have to be regularized in the UV (e.g. dimensionally) 
and renormalized such that, e.g., the renormalization conditions 

/i,a(0) = 1 , / 2 ,a(0) = 1 (4.3) 
are satisfied. As a result one obtains 

/i,a(p 2 ) = 1 + 5/ 1 1 - 1o ° p (p 2 ) , 

/ 2 ,a(p 2 ) = 1 + ^ 1o ° P (p 2 ) (4.4) 
where the one-loop contributions 5fl~ loop are given in the appendix. For 
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large p 2 the functions 5fl loop behave asymptotically as 
6 - 




5/-V)~(^-iln^j ■ (4.5) 

What we have achieved now is an appropriate behaviour of these functions 
even for momenta p 2 3> A 2 , far above the starting scale, provided g\j (47r) 2 ■ 
ln(p 2 /A 2 ) <C 1. Now Vk(p 2 ) of eq. (3.14) has the appropriate logarithmic in- 
dependence even for p 2 ^> A 2 , corresponding to a one-loop RG improvement 
of the "physical" coupling constant: 

\4 =A (p 2 »A 2 ) -4 . (4.6) 




As boundary conditions for the integration of the ERGEs we thus use 
propagator functions f it \ according to eqs. (4.4), the numerical solution of 
eq. (3.15) for m\, and values for (as the only free parameter) from g\ = 2.0 
down to (?a = 1-4. 

It is clear that, in order to integrate the ERGEs (4.1) for the functions 
fi t k and /2,/fc, numerical methods have to be employed. To this end we need 
a parametrization of these functions, which allows for a reliable fit for all 
values of the scale k. Since these functions always behave like fi,k{p 2 
oo) — > 1 +5/j l loop , the following parametrization turns out to be convenient: 

/^^l-E-b + r° P (/) (4.7) 

j h r 

with up to 6 pairs of parameters a*, 7* for % — 1, 2, respectively. 

Next we have to diagonalize eqs. (4.1) in the space of momenta p 2 . To 
this end we evaluate the functions h\ in eqs. (4.1) and (4.2) for a number of 
momenta p 2 = p 2 with up to 35 different values for p 2 in the range lO -1 /^ 
to 10A 2 , where /cq is determined dynamically (see below). (It turns out to 
be convenient to space the momenta p 2 equally on a logarithmic scale.) By 
solving the system of linear algebraic equations we then obtain expressions for 
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dkfi,kipf) an d d k m 2 k . One step in the integration of the ERGE then amounts 
to a computation of 

/i,k+Afc(p?) = fi,k{pj) + AA; • d k fi, k (pj) (4.8) 

with \Ak/k\ < 3 • 1CT 2 , and a subsequent determination of the new fit pa- 
rameters a*-, 7] of eq. (4.7). 

The new value of the parameter m\ gets determined by a new numerical 
solution of eq. (3.15) rather than by the numerical integration of d k m\\ since 
it is a relevant parameter (the only one), for k 2 <C A 2 accumulated theoretical 
and numerical errors would easily generate a m\ of 0(A 2 ) instead of 0(k 2 ), 
as it should be according to eq. (3.15). 

Of course we have varied many details of the numerical procedure in order 
to test the robustness of the final results. An upper limit of ~ 35 different 
values for p 2 is dictated by the finite amount of available computing time, 
and this also keeps us from choosing arbitrarily small values for the gauge 
coupling g A at the starting point k 2 = A 2 , since then the ERGEs would have 
to be integrated over arbitrarily many orders of magnitude of the scale k 2 in 
order to arrive at the physically interesting non-perturbative regime. 

5. Results 

Our results can most easily be represented in terms of a function F k (p 2 ), 

F k (p 2 ) = ± fiAP 2 ) fl(P 2 ) , (5-1) 
9a 

which is simply related to the potential V k {p 2 ) in eq. (3.14): 

V k (p 2 )= 1 • (5.2) 
p 2 F k (p 2 ) 

In figs. 2-4 we show our results for F k (p 2 ) for a bare coupling g\ = 1.4. We 
plotted this function for different values of k 2 , k 2 = A 2 , 10 _1 A 2 , 10~ 2 A 2 , 
1(T 3 A 2 , k 2 = k 2 = 7.43 ■ 1(T 5 A 2 , as full lines, and for some still smaller 
value k 2 = k = 3.68 ■ 1CT 5 A 2 as a dashed line, p 2 is given in units of A 2 ; 
the maximal value of p 2 is A 2 , 10~ 1 A 2 and, in order to resolve the small 
p 2 region, 10 _2 A 2 in figs. 2 to 4, respectively. The curves with decreasing 
values at p 2 = correspond to decreasing values of k 2 . We see that, for 
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small enough k 2 and small enough p 2 , F k (p 2 ) approaches a form ~ p 2 near 
the origin, which corresponds to a l/p 4 -behaviour of the potential V k {p 2 ). 

Just before F k (p 2 ) becomes for p 2 — > 0, however, the two different 
methods of evaluating d k m\ (cf. eq. (4.2) and the discussion below) lead to 
different results for d k fi tk or d k f2 tk , with a relative difference of ~ 0.1. We 
denoted the corresponding scale by fc . For slightly smaller values of k the 
equation for the gluon mass term m 2 k (3.15) possesses no longer a solution 
for m 2 k . 

The incompatibility of the RG flow with the modified STIs (3.8) be- 
low k indicates, that our approximation becomes unreliable in this regime: 
Numerically different results for d k f\ jk or d k f2,k show, that the neglected con- 
tributions play an important role on the r.h.side of the ERGE for the gluon 
or ghost propagator, respectively, and that our ansatz is too restrictive in 
this regime. 

The gluon mass term m 2 k itself, which was determined by eq. (3.15) to be 
of 0(A 2 ) at the starting point k = A, does not show an unusual behaviour 
for k ~ k . ml is always of 0{k 2 ) or, more precisely, m\ ~ —0.6 • k 2 for 
k — > k . 

Thus we cannot have confidence in the form of the RG flow for k < 
ho. However, generally the RG flow modifies F k (p 2 ) only for p 2 < k 2 (cf. 
figs. 2-4), hence our result for F ko is nevertheless trustworthy for p 2 > ko . 
Our results also indicate that, within a more general parametrization of f k , 
F k (0) approaches for k 2 — > 0: In figs. 2-4 we have shown the function 
F k (p 2 ) for some value k 2 = k 2 < k 2 as a dashed line, and we see that the 
decrease of F k (p 2 ) at p 2 = continues. For still smaller values of k 2 , however, 
the integration of the system of differential equations becomes numerically 
unstable and does no longer allow to obtain reliable results. 

In any case, already our results for F ko (p 2 ) for p 2 > k 2 strongly indicate 
a form ~ p 2 near the origin, and the trustworthy range of p 2 may be sufficient 
for some phenomenological investigations involving, e.g., heavy quarks. 

It might be interesting to get some feeling for the orders of magnitude of 
the different scales involved. Up to now all scales are only defined relative to 
the starting scale A. In order to relate this scale to a physical scale one can 
try to compare our result for V ko (p 2 ), for p 2 > k 2 , with a phenomenological 
parametrization of this potential, which contains dimensionful parameters 
known in units of MeV. A convenient form of such a parametrization is given 
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by Richardson [24], which reads in our convention 



Vr( P 2 ) - (33 _ 2Nf)p2 ln (1 + p2/A2R) (5.3) 

with Ar = 400 MeV, and for a rough comparison with our results (in pure 
SU(3) Yang-Mills theory) we may set Nf = 0. In terms of our function 
F k (p 2 ) of eq. (5.1) this ansatz reads 

= ilMigVAl) . ( , 4) 

We performed a fit of eq. (5.4) to our function F ko {p 2 ) for the region of 
momenta k 2 , <C p 2 <C A 2 to determine the ratio A/A R . The result is shown 
in figs. 5 and 6, where we plotted F ko {p 2 ) as a full line and FrQo 2 ), eq. (5.4), 
as a dashed line, against p 2 in GeV 2 . The optimized relation between A and 
Ar is such that A = 12.1 GeV, which implies that the scale ho, where our 
approximation ceases to be appropriate, is ~ 104 MeV. One can see that our 
potential is surprisingly close to the phenomenological one of Richardson, 
both in the perturbative (fig. 5) and in the nonperturbative (fig. 6) regime. 
Significant deviations show up only for momenta p 2 < . 

Let us now look at the dependence of our results on the bare coupling g^. 
We find that also for larger values of g\, F k (0) becomes small for small k, until 
our consistency condition indicates the inadequacy of the approximation at 
some scale k , with k 2 /A 2 = 6.13-10" 4 , 2.69-10~ 3 , and 7.98-10' 3 for g A = 1.6, 
1.8 and 2.0, respectively. In fig. 7 we compare the different results for Fk (p 2 ) 
with g\ = 1.4 (full line), g\ = 1.6 (long-dashed line), g\ = 1.8 (short-dashed 
line) and g\ = 2.0 (dotted line). In each case we determined the physical 
value of the starting scale by a fit to eq. (5.4), and the results for A are 
given by 12.1 GeV, 4.8 GeV, 2.6 GeV and 1.7 GeV for g A = 1.4, 1.6, 1.8 and 
2.0, respectively. From the values of k^/A 2 given above, we then obtain the 
corresponding values for k , k = 104 MeV, 119 MeV, 134 MeV and 151 MeV. 
We see that the curves in fig. 7 coincide remarkably well, which assures us 
that the numerical method we use for the integration of the ERGEs does not 
accumulate numerical errors. The approximate independence of the curves 
on g\ or A even for momenta larger than the respective starting scale is due 
to the one-loop improvement of the starting action T\ (or of the functions 
/i 5 a and / 2; a according to eq. (4.4)). In fig. 8 we show the same curves for 
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Fk {p 2 ) f° r the different values of g\ or A, which have been obtained without 
the one-loop improvement. 

Figs. 7 and 8 lead to a series of important observations: First, as we 
see most clearly from fig. 8, the curves are essentially unchanged compared 
to their form at the respective starting scale A for momenta p 2 ^> A 2 , as 
they should be. In particular, the curve for A = 1.7 GeV (g\ = 2.0) in 
fig. 7 is given in the perturbative region essentially by the one-loop form 
(4.4), so comparison with, e.g., the curves for A = 12.1 GeV (g\ = 1.4) in 
both figures shows, that the integration of the ERGEs correctly reproduces 
the one-loop form in the perturbative domain. More importantly, we find 
universal behaviour in the non-perturbative domain: For momenta p 2 < A 2 , 
i.e. below the respective starting scales, all curves merge in one, as is most 
impressively demonstrated by fig. 8, and the form of the resulting curve does 
not depend on the boundary condition we use. This independence of physics 
at p 2 <C A 2 with respect to variations of the starting scale A and to the 
inclusion of irrelevant operators at A (as is the case for the one-loop improved 
boundary conditions) is what we call universality here, in accordance with 
ref. [23]. 

It is easy to see that the relation between and the starting scale A, 
which has been obtained empirically above, agrees with the one-loop running 
of g\. On the other hand the running of g\ with A is not yet correct to two- 
loop accuracy. However, with our ansatz for only a subset of all possible 
two-loop diagrams is included within a perturbative expansion of the ERGEs, 
hence an agreement at this level would have been mere luck. 

One may ask whether the behaviour of Fk (p 2 ) for p 2 — > is due to the 
form of fi t k {p 2 ) or both fi,k Q (p 2 ) and f2,k {p 2 )- To this end we plot in fig. 9 
both functions fi,k (p 2 ) an d f2,k (p 2 ) (with p 2 rescaled as before) for g\ = 
1.4. We see that both functions become small for p 2 — > 0, but it would be 
grossly misleading (and unphysical) to identify the gluon propagator function 
(p 2 /i(p 2 )) -1 with the heavy quark potential and to neglect the dependence 
of V(p 2 ) of eq. (3.14) on f2{p 2 ), as it is done if the ghost part of the action 
is neglected, and the STIs are not properly taken into account. 

6. Discussion and Outlook 

In this paper we have performed the integration of the ERGEs for non- 
abelian gauge theories, with the bare action as the only input. We con- 
centrated on pure SU(3) Yang-Mills theory (or quenched QCD), and the 
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momentum dependence of the gluon and ghost propagators; these quantities 
determine the heavy quark potential. 

The neglect of certain contributions on the r.h. sides of the ERGEs for 
the 2-point functions turned the system of ERGEs into a closed system of 
differential equations, which we integrated numerically. A priori the validity 
of the approximation is difficult to judge in the strong coupling regime, where 
it amounts to the neglect of contributions of the same order in the coupling 
constant as the ones maintained. Here the use of modified STIs plays an 
important role: Due to their validity parts of the /c-dependent effective ac- 
tion can be determined in two different ways, namely either directly through 
these identities or by integrating the ERGEs. Without approximations both 
methods give the same results, but in the presence of approximations a pos- 
sible difference between the corresponding results can be used to estimate 
the error induced by the approximations. In our case it was the gluon mass 
term, which was obtained in two different ways. We used a 10%-deviation 
in dkfi,k or dkf'2,k, induced by the different methods to obtain dkm 2 ., as an 
indication for the breakdown of our approximation near the corresponding 
scale k . Accordingly, our result for V(p 2 ) = Vk (p 2 ) is only consistent for 
P 2 ~ k 2 , where we estimated k ~ 100 MeV. Our main findings are the 
strong indication of a 1/p -behaviour of V(p 2 ) in the trustworthy regime of 
p 2 , and the fact that the full form of V{p 2 ) happens to be close to the one 
proposed by Richardson [24]. 

In the second section of this paper we discussed the formal relation be- 
tween ERGEs and SDEs, with the result that effective actions, which solve 
SDEs in the presence of an IR cutoff k, can be considered as quasi-fixed 
points of the ERGEs. Although this formal relation is no longer exact in the 
presence of truncations of effective actions, it is thus not astonishing, that 
the behaviour of the gluon propagator, which has been derived here using the 
ERGEs, has also been shown to be a solution to SDEs in the Landau gauge 
[17]. We believe, however, that the present method possesses a number of 
advantages compared to the formalism of the SDEs: Firstly, the regulariza- 
tion of both ultraviolet and infrared singularities in the context of the SDEs 
generates notorious technical difficulties, which require, to some extent, ad 
hoc prescriptions, the influence of which on the final result is difficult to con- 
trol. The present method is, in contrast, free of problems related to UV or 
IR singularities by construction. 

Secondly, a priori both methods require a suitable truncation of the ef- 
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fective action. Within the present method, however, it is technically simpler 
to include more and more terms; for the present investigation, e.g., we al- 
ready included the contributions of the 4-gluon vertex and the ghost fields, 
which were neglected in the SDE approach [17]. We have stressed that the 
detailed form of the dressed one-gluon exchange diagram between heavy 
quarks, which involves the ghost propagator function at the quark-gluon 
vertex due to the STIs, differs considerably from the form of the inverse 
gluon propagator alone (in the Landau gauge), which is not a physical quan- 
tity by itself. Moreover, it is technically feasible in our formalism to include 
further terms in the effective action without too much effort. 

Finally, within the SDE approach one is only able to check the self- 
consistency of a particular ansatz; it is practically impossible to prove rigor- 
ously, that the chosen ansatz is the only possible one. The integration of the 
ERGEs, on the other hand, always gives a unique answer, namely the one 
which is smoothly connected to the corresponding bare action, and to per- 
turbation theory in the UV in the case of asymptotic freedom. Actually, our 
discussion in section 2 sheds some light on the cases where different solutions 
of the SDEs exist: on the one hand solutions of the SDEs can be considered 
as fixed points of the ERGEs for k — > 0, on the other hand it is generally not 
clear, whether one has obtained an IR stable or an UV stable fixed point. 
Since only the IR stable solutions can also be obtained by the integration of 
the ERGEs towards k — > 0, one is lead to the conclusion, that only these are 
actually physically relevant, which is certainly difficult to check within the 
SDE approach alone. 

Many possibilities exist in order to extend the present approach in the 
future: firstly, still more terms can be included in the effective action, and 
collective fields, e.g. for F® U F® U , can be introduced along the lines of refs. [7, 8] 
in order to reach renormalization scales k below ho, and to get a still better 
understanding of pure Yang-Mills theory free of systematic uncertainties 
for p 2 < /c 2 . Secondly, the quark sector of QCD can be introduced as 
well. First efforts in this direction have already been made in ref. [8]; there, 
however, the effects of the gluons have been guessed and parametrized in a 
phenomenological manner. Now the integration of the ERGEs for the full 
quark gluon system is within reach (possibly along the lines proposed in ref. 
[25]), which will allow to study many phenomenologically interesting systems 
like heavy gg-bound states, gg-condensates and the meson sector on the basis 
of just the bare QCD Lagrangian. 
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Appendix 

In this appendix we list some of the functions which are needed in the 
actual computations. We begin with the 3- and 4-gluon vertices employed 
in our ansatz (3.10) for rV They are chosen in such a way as to fulfill the 
standard STIs in a "minimal" fashion. Explicitly they read 



and 



(qjg,k(<l 2 ) - r^fgAr 2 )) K + (r u f g Ar 2 )-Puf 9 Ap 2 )) 

+ (ppfgAp 2 ) - QpfgAq 2 )) V 

+ {% - rp) UAq2 yj^ {r2) M, - q ■ r 5 vp ) 

+ (r„ - p.) f9 ^lZp2 k{p2) (Pp^ -r-P V) 
+ (P P - q P ) f ^ p2 yj^ q2) { %Vv -p.q 5 pu ) | (A. 1) 



W(p,9,r,s) 4/ 2 , fe (^) /2;fc ( ? 2 )/2ifc(r 2 )/2ifc(s2) >< 



/p,*((p + ?) 2 )^p <W - 2 - ^^^^ ( P(7 S p -p. 8 5 up 

p s 



(r + s) 2 — s 2 \ p 2 — (r + s) 2 p 2 — s 2 J 

x (2p u + q u )(2s p + r p )(p IT s p -p- s5 ap )^j , (A.2) 

where 

W) = /i, fc (p 2 )/ 2 2 , fc (p 2 ) • (A.3) 
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Observe that we have not symmetrized the expression (A. 2) for the four- 
gluon vertex with respect to momenta and indices, in order to keep it to a 
reasonable length. 

Next we turn to the modified STI for the gluon mass term, eq. (3.15). In 
diagrammatic form it is given by a number of one-loop diagrams. The inte- 
grations over the angular variables of the loop momenta can be performed, 
and one is left with integrals over t = q 2 . After a careful consideration of the 
limit a — > 0, corresponding to the Landau gauge, the following expression 
for the function ST results: 



ST '(k 2 , ml, /i,fc,/ 2) fc) 



4/1. 



oo ✓ 



t d 



+ ( 3 -^m( thAt) )) f2At) 



1 _ e -{th,k{t)+m 2 k )/k' 2 

tfi,k(t)+m 2 k 

_ p -tf2,k(t)/* 2 ~ 



tf2,k(t) 



_l_ 3 e -(*/i,fcW+™l)/fc 2 



' 4 -4|(*/m(*))V i .*w- 



-(tf 1>k (t)+m 2 )/k 2 



tfi,k(t) +m 



+ f2,k(t) 



1 _ e -t/ 2 , fc (t)/fc 2 
t/ 2 ,fc(t) 



6 RT^ x 



/ 2lfc (*) <9t 



(t/2,*(* 



1 _ e -(tfi,k(t)+m 2 k )/k 2 

tfi,k(t)+ml 



_ t d_( tfik{t) j 1 ~ (l + + m 2 )/k 2 )e-^)+rn 2 )/ k 2 



(tf 1>k (t)+ml 



(A.4) 



Finally we give the one-loop expressions 5fl~ loop for the gluon and ghost 
2-point functions as needed in (4.4) for the improved boundary conditions. 
They are calculated with the IR cutoff functions R k from eqs. (3.16) (with 
and f 2; k set equal to one and m\ = 0) and obey the renormalization 
conditions (4.3). The results are, again in the Landau gauge, 



= H{lf(M*/2) + c)-H + 



27x 2 + 88x - 20 
24^3 



28 



x 4 - 8a; 3 + 9x 2 - 20s - 6 



and 



12a; 3 



e x + 



x 4 — 16x 3 — 16x 2 — 56x + 4 
12x3 



x 2 - 7x _. . , a; 2 - 14a; - 52 ^. . . . ' 
Ei(-x) H — Ei(-x/2) 



12 



24 



rr — Z- 

A 2 



5ft loop ( P 2 ) 



3gj 
(4tt) 2 



|(ln(,/2) + C)-| 



3x + 2 x 2 - x + 2 
+ : o + e 



x 2 + x + 2 



-x/2 



4x 2 



a; 



4x 2 
x + 3 



2x 2 



+ I Ei(-x)^ ( 



Ei(-x/2) 



A 2 



-x/2 



(A.5) 



(A.6) 



where C denotes Euler's constant, C = 0.577216. The exponential integral 
function Ei is defined by 



~ x t 
Ei(-x) = J dtj 



(A.7) 



for x > 0. 
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Figure captions 

Figure 1: Diagrammatic form of the ERGEs for the ghost propagator func- 
tion / 2) fc and the gluon propagator function /i^. Internal curly lines 
denote the full gluon propagator, internal dotted lines the full ghost 
propagator and full points the full vertices from the ansatz (3.10) for 
T k . The crossed circles denote insertions of d^Rk^u resp. d^Rk- 

Figure 2: Results for the k-dependent function Fk(p 2 ) (5.1), which is re- 
lated to the potential Vk(p 2 ) through (5.2), for different values of k 2 , 
k 2 = A 2 , lO^A 2 , 10~ 2 A 2 , 10" 3 A 2 , k 2 = k 2 = 7.43-10" 5 A 2 , as full lines, 
and for k 2 = k = 3.68 • 10~ 5 A 2 as a dashed line. The bare coupling g A 
is (?a = 1.4, and p 2 is given in units of A 2 . The curves with decreasing 
values at p 2 = correspond to decreasing values of k 2 . 

Figure 3: As in fig. 2, with a maximal value of p 2 = 10^ 1 A 2 . 

Figure 4: As in fig. 2, with a maximal value of p 2 = 10 _2 A 2 in order to 
resolve the small p 2 region. 
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Figure 5: Fit of our result for Fk (p 2 ) (full line) to the parametrization of 
Richardson [24] (dashed line), with an optimized relation between A 
and A R ^ 400 MeV such that A = 12.1 GeV. 

Figure 6: As in fig. 5, with a maximal value of p 2 = 4 GeV 2 . 

Figure 7: Results for Fk (p 2 ) for different values of the bare coupling g\ = 
1.4 (full line), g\ = 1.6 (long-dashed line), g\ = 1.8 (short-dashed 
line) and g\ = 2.0 (dotted line). In each case we determined the 
physical value of the starting scale A by a fit to the parametrization of 
Richardson, eq. (5.4). 

Figure 8: As in fig. 7, but without the one-loop improvement of the prop- 
agator functions at the starting scale A. 

Figure 9: Results for the propagator functions f\ t k (full line) and $2,k 
(dashed line) for g\ = 1.4. 
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